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Abstract. Identities and inequalities are proved for the order parameters, correlation 
functions and their derivatives of the Ising spin glass. The results serve as additional 
evidence that the ferromagnetic phase is composed of two regions, one with strong 
ferromagnetic ordering and the other with the effects of disorder dominant. The 
Nishimori line marks a crossover between these two regions. 
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1. Introduction 

One of the outstanding problems in the theory of finite-dimensional spin glasses is the 
structure of the phase diagram. Numerical investigations have revealed the precise 
locations of critical points and phase boundaries However we still have only limited 
knowledge from analytical treatments of the problem. An interesting exception is the 
gauge theory which makes use of gauge symmetry of the system to derive a variety 
of exact /rigorous results on physical quantities including the energy, specific heat and 
correlation functions || |3|. In the present paper we derive a class of relations for 
the order parameters and correlation functions using the gauge theory to clarify the 
behaviour of the system within the ferromagnetic phase. 

Properties of the ferromagnetic phase in models of spin glasses have not been 
studied very extensively compared to the spin glass phase. Nevertheless, as shown 
in references [|J and ||, a very interesting non-trivial change of the system behaviour 
is observed on a line, the Nishimori line, in the phase diagram: The spins become more 
ferromagnetically ordered (i.e. parallel to each other) as the temperature is lowered 
above this line whereas the same spins turn to become misaligned below the same line 
with further decrease of the temperature. Although this line is not a phase boundary 
in the thermodynamic sense, it marks in the above sense a clear crossover between the 
two regions within the ferromagnetic phase. The argument in the present paper using 
the gauge theory reinforces this picture through relations between the order parameters 
and their derivatives. 

An important consequence of the gauge theory is an identity between the 
distribution function of the ferromagnetic order parameter P m (x) and that of the spin 
glass order parameter P q (x): These two functions are equal to each other P m (x) = P q (x) 
on the Nishimori line |3], |J. This functional identity implies the absence of replica 
symmetry breaking because P m (x) is trivial, composed of at most two simple delta 
functions, and, therefore, so is P q (x). The relation P m (x) = P q (x) also leads to the 
equality m = q, an identity between the ferromagnetic and spin glass order parameters 
0, H , implying an exact balance between the two types of ordering. From m = q we may 
expect that the ferromagnetic order parameter exceeds the spin glass order parameter 
m > q above the Nishimori line because, in the limit of non-random system (which 
is above the line), we have q = m 2 < m. Another reason to expect m > q is that, 
as mentioned above, ferromagnetic ordering dominates above the line. The opposite 
inequality q > m is likely to hold below the same line since the spin glass phase (lying 
below the line) has q > m = 0, and, in addition, the effects of quenched randomness is 
more dominant (apparent misalignment of spins) below the line as explained above. 

These two inequalities for the order parameters can be verified in the mean-field 
Sherrington-Kirkpatrick model J7| within the replica-symmetric solution (which is valid 
near the Nishimori line as mentioned above) but have been considered difficult to check 
analytically for general finite-dimensional systems. We show in the present paper that 
these inequalities are closely related with temperature derivatives and some inequalities 
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Figure 1. A plausible generic phase diagram of the ±J model. The system has 
paramagnetic (P), ferromagnetic (F) and spin glass (SG) phases. The Nishimori line 
K = K p shown dashed marks a crossover between two regions within the ferromagnetic 
phase. 



of the order parameters, thus coming closer to a proof. 

We present our results and their proofs in the next section. Discussions are given 
in the last section. Some of the details of calculations are described in the Appendix. 

2. Identities and inequalities 

To be specific, let us consider the ±J Ising model on an arbitrary lattice with the 
probability of ferromagnetic interaction denoted by p. The expected phase diagram 
is depicted in figure [p. The main physical quantities we treat in this paper are the 
ferromagnetic and spin glass order parameters defined by 

m = [(S i ) K ], q=[(S i ) 2 K ], (1) 

respectively, with i well interior of the lattice. The inner brackets denote the thermal 
average with coupling constant K = J/T, and the outer square brackets are for the 
configurational average characterized by the parameter p. The spins on the boundary 
of the lattice under consideration are set to up states to avoid trivial vanishing of the 
thermal expectation value of the one-point function (Si)k- It is sufficient to consider a 
single parameter q as a spin glass order parameter instead of its distribution function 
P q (x) because replica symmetry breaking is absent when K = K p = |log(p/(l — p)) 
H |3|], the Nishimori line, on and near which we concentrate ourselves for the moment. 

The identity m = q valid when K = K p has long been known [[|. The first of our 
new results is the following relations: 
dq_ = dm 

dK dK 1 ' 

OK 2 ~ dK 2 ' { ' 

both of which hold under the condition K = K p . 
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The proof is straightforward. As shown in the Appendix, the magnetization is 
rewritten using gauge transformation as 

m=[(S i ) K ] = [(cr i ) Kp {S i ) K ], (4) 

where (o~i)K p is the thermal average of the Ising spin (introduced by the gauge 
transformation) of the same system as the original ±J model with effective coupling 
K p . The first and second derivatives of m are then 



dm 
~dK ~~ 
d 2 m 
dK 2 



d 

\Vi)K p -Q^{Si)K 

d 2 



A" 



(5) 



(6) 



The derivatives of the spin glass order parameter are obtained directly from the definition 



dq 

~dK ' 

d 2 q 
dK 2 



\Si)K 
{Si)K 







dK 

Q2 



dK 2 



{Si)k 

(3i) K 



+ 2 



' d_ 
dK 



(Si) 



K 



(7) 



(8) 



The identity m = q for K = K p immediately follows from and (|T|) because 
(vi)Kp = (Si)x if K = Kp. The identity (||) valid for K = K p is also easy to verify from 
(|) and (|7]). The inequality (||) is a consequence of (|]) and (|^). 

Similar relations hold for more general correlation functions. Let us define two 
correlations: 



M21+1) _ r/Q Q C \2/+l] M21+2) _r/q q q \2l+2] 



where / is a non-negative integer and {i,j, k, ■ ■ •} is an arbitrary set of sites. 

y(2Z+l) 



This C^Ti^y satisfies the following identity (see the Appendix) 



c ifk + - 1] = [(° r i°~j°~k ■ ■ ■)K p (S l S j S k ■ ■ ■) 



21+1] 
K J 



(9) 



(10) 



Using the fact that C\f^ is gauge invariant, we can prove the following relations at 
K 



Kp. 



n {2l+l) 
^ijk- 



M21+2) 
^ijk--- i 



1 9 r (2l+l) 

21 + 1 OK ijk - 



1 



d 



M21+2) 

21 - 2 dK ijk - ' 



No simple relation exists between the second derivatives for general I. It is to be noted 
that equation ( |Tl"D holds not just in the ferromagnetic phase but in the paramagnetic 
phase as well whereas equations (0) and (^|) are trivial in the paramagnetic phase as 
both sides vanish identically. 

We next discuss our second new result for the order parameters. Let us denote 
the dependence of the order parameters on the temperature and probability parameter 
explicitly as m(K, K p ) and q(K, K p ). Then it is possible to show that 

m(K,K p ) >q(K,K p ) 

=>- m{K p , K p ) > m{K, K p ) and q{K p , K p ) > q{K, K p ) (12) 
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for any K and K p . To prove this, it is useful to write (0) explicitly as 



m{K,K p ) =E^(*)fo>g<S>K = £ • y/PWiSiW, (13) 

where stands for a bond configuration. Let us square both sides of the above equation 
and apply the Schwarz inequality to obtain 

" ' \2 



i(k,k p ) 2 < e^)(wS) 2 E^)((^)? ) : 



m( , _ 

A- 

v2 ir/o\2 

q(K p ,K p )q(K,K p ). (14) 



Now, if we assume m(iT, Jf p ) > iT p ), then g(if, K p ) on the right hand side can be 
replaced by m(K, K p ) to yield 

m(K, K p ) 2 < q(K p , K p )m(K, K p ). (15) 

Since we are considering the ferromagnetic phase with up-spin boundaries, we have 
m(K, K p ) > 0, and therefore by dividing both sides by m(K, K p ), we find 

m(K, K p ) < q{K p , K p ) = m{K p , K p ). (16) 

This is the first half of the result (P^). 



The second half is proved similarly. From fllj) and the assumption m(K } K p ) > 
q(K, K p ), we find 

q(K, Kp) 2 < q(K p , K p )q(K, K p ) (17) 

and thus 

q{K,K p )<q(K p ,K p ). (18) 
A generalization to correlation functions is straightforward. The result is 

C\l + }\K,K p )>C^ +2 \K,K p ) 

=► Cg..(K P ,K p ) > C$£!>{K,K P ) and Cg..(K p ,K p ) > c\% +2 \k,K p ). 

(19) 

To prove these inequalities, we apply the gauge transformation and Schwarz inequality 
to C\% + }\K,Kp): 



Cg + }\K,Kp) 2 =[{S l S J S k ...) 2 ^} 2 



= (E y/PWiWjVk ■ ■ ->Sg ■ ({StSjSk ■ ■ -)P) 2l+ ^ 

k k 

= [(o-iO-jO-k ■ ■ ■) 2 Kp ][(S l S J S k ■ • -) 4 K +2 } 

= C%UKp,K p )C% l k +2 \K,K p ). (20) 
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Figure 2. It is forbidden that the derivatives of m and q have different signs at 
K = K p as shown in this figure. 



When C\f£\K, K p ) > Qj k Z } (K, K p ), the second factor on the right hand side of (0) 
is bounded from above by C^ l k + } (K, K p ) to yield 



the final inequality being a result of gauge transformation of the kind described in the 
Appendix. This ends the proof of the first inequality of (0). By replacing the left hand 
side of © with the lower bound C^ 2 \k, K p ) 2 , we arrive at the second relation 



It is possible to apply similar arguments to the other models of spin glasses 
with gauge symmetry including the random Ising model with Gaussian-distributed 
interactions and XY gauge glass || |J . The physical significance of the results obtained 
in this section will be discussed in the next section. 

3. Discussions 

An immediate consequence of @ is that the derivatives of q and m have the same sign 
when K = K p . It is forbidden that, for example, the spin glass order parameter q 
increases whereas the ferromagnetic order parameter m decreases as depicted in figure 
H]. In the plausible case that dq/dK > 0, it follows from @ that q increases twice 
as rapidly as m does as the temperature is lowered (see figure |3| (left)). It naturally 
follows that q is larger than m at least slightly below the Nishimori line (K > K p ) and 
the opposite inequality m > q holds when K < K p . This is the most natural 
discussed in the introduction: The ferromagnetic ordering is dominant (m > q) above 
the Nishimori line (K < K p ) and quenched-disorder-driven random ordering proliferates 
in the sense q > m below the same line (K > K p ). The point is that we have reduced the 
inequality q > m for K > K p to an intuitively natural relation dq/dK > at K = K p 
(although we still do not have a rigorous proof of the latter.) 




(21) 



(22) 



C\%Z Z \K ) K p )<Clf k ...{K p ,K p ). 



(23) 
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Figure 3. The spin glass order parameter increases twice as rapidly as the 
ferromagnetic order parameter around K = K. p if these quantities are increasing 
functions of K (left). The opposite possibility is shown on the right. 



If the opposite inequality dq/dK < holds at K = K p , the order parameters 
behave as depicted in figure |3] (right). We cannot exclude this possibility from the present 
argument, but it seems quite unlikely that both order parameters q and m decrease with 
temperature decrease on and around the line K = K p which runs through the central 
part of the ferromagnetic phase as seen in figure |l|. If there exists a reentrant transition 
(below which m vanishes but q can stay finite), only m may decrease toward such a 
transition temperature from above unlike figure (right). Although we believe that such 
a situation is not plausible to exist 0, it would happen at very low temperatures if it 
does at all, not around the line K = K p which runs through relatively high temperature 
parts of the ferromagnetic phase. 

The final possibility is that the derivatives in (Q) vanish. Again, both derivatives 
should vanish, not just one of them. The inequality for the second derivative (0) does 
not tell much about the behaviour of m and q around K = K p . 

It is instructive to remember in this connection that the average sign of local 
magnetization reaches its maximum at K = K p as a function of K (or the temperature) 



M(K, K p ) = [sgn(^] < M(K p , K p ). (24) 

This means that the number of up spins under up-spin boundaries becomes maximum at 
K = K p as a function of K. Although this relation (|2~4f) can be proved for an arbitrary 
K by the gauge theory || , it is useful to check it by taking the derivative of 
did for m and q): 



as we 



dM 
~dK 



_d_ 

dK 



sgn(Si) 



A' 



K) 



d 
~dK 



(Si) 



K 



K , 



d 
~dK 



(Si) 



K 



0, (25) 



the last equality being valid for K = K p . Equation (p5|) implies that the effects of spins 



with positive temperature derivative (d(Si)K/dK > 0) just match those of negative 
temperature derivative (d(Si)K/dK < 0) at K = K p if we concentrate ourselves on 
the spins with vanishing local magnetization S^S^k)- Thus, at K = K p , some spins 
change its local magnetization from positive value to negative value whereas essentially 
the same number of spins change the sign of local magnetization in the opposite way. 
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It should be noted, however, that this observation does not necessarily mean vanishing 
derivatives of the ferromagnetic and spin glass order parameters at K = K p , that is, 
dm/ OK and dq/dK are in general not vanishing at K = K p . The absolute value of 
local magnetization | which is ignored in M(K), at sites with > may 

grow more rapidly than at sites with (Si)x < as K decreases, compensating for the 
decrease in the number of up spins below K = K p , leading to a positive derivative 
dm/dK > 0. 

The present argument does not apply in the paramagnetic phase where the order 
parameters m and q vanish. However, the relation (|TT|) for correlation functions, 
the two-point functions C^ l+r> (K, K p ) and C^ l+2 ' (K, K p ) in particular, suggests that 
C\f +1 \K,K P ) > C\f +2 \K,K P ) if K < K p and C\f +l \K,K p ) < Cf +2) {K,K p ) 
if K > K p . This means that the ferromagnetic correlation length £ m (defined by 
C$ ~ exp(— r/£ m ), r ^> 1) is larger than the spin glass correlation length £ q (defined by 
CqP ~ exp(— r/£ ? ),r ^> 1) for K < K p (above the Nishimori line) whereas the opposite 
inequality holds below the line. 

Very similar conclusions follow from the result (0). The relation between m 
and q shown in figure |2| violates these inequalities because, in the temperature range 
where m(K,K p ) exceeds q(K,K p ), m(K,K p ) is seen to be larger than its value at 
K = K p . The cases given in figure |3| are compatible with the present inequalities: If 
m(K, K p ) > q(K, K p ), then m(K, K p ) is smaller than m(K p , K p ). An advantage of the 
present inequalities ([12]) over the relation (Q) is that we can analyze the behaviour of 
order parameters well away from the Nishimori line, that is, K can be much larger or 
smaller than K p . A weak point is that quantitative relations for the increase/decrease 
rates are not given unlike (||). Another problem to remember concerning ([12]) is that 
we may not be able to describe the system only in terms of m and q at very low 
temperatures if replica symmetry breaking exists as it is the case in the Sherrington- 
Kirkpartrick model [TO . 

The analysis presented above strongly indicates that the Nishimori line K = K p 
marks a crossover between the purely ferromagnetically-ordered region and disorder- 
dominated region within the ferromagnetic phase. This observation is consistent also 
with renormalization group analyses: The low-temperature region is controlled by a 
fixed point at T = with finite disorder (p < 1) whereas the high-temperature region 
is described by a different fixed point representative of the critical curve above the 
multicritical point [11, |T^] . 

An important future direction of investigation is a rigorous proof of the relation 
dm/dK > at K = K p , which needs additional ideas. 
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Appendix 

In this Appendix we drive equations (|) and fllUD following references |2| and ||. The 
ferromagnetic order parameter is defined by 

~ [{Sl)Kl ~ (2cosh K 9 )»* {T ^ ±l} M P h l3) E { 5 } exp(AE fo> r^) " (26) 

Here t^ is the sign of J y - (r^ = Jij/\Jij\ = ±1), the sums in the exponents run over all 
pairs of sites on the lattice under consideration, and N B is the total number of bonds. 
The factor exp (K p Tij)/2 cosh K p gives the probability weight of configurational average 
of the ±J model since this quantity equals p for Ty = 1 and 1 — p for Ty = — 1 as can 
be verified from the definition K p = | log(j»/ (1 — p)). Spins on the boundary are all up. 
Let us apply a gauge transformation 

S{ > Si&i, Tij ► TijO~iOj (27) 

to all sites, where o~i is a gauge variable fixed either to 1 or —1 arbitrarily at each site 
(+1 on the boundary). The Hamiltonian in the exponents of fl26|) is invariant under this 
gauge transformation. Since the gauge transformation is just a re-definition of running 
variables in (^), it does not affect the value of the right hand side, and we have 

(2coshA p )^ { ^ } eXp ^^ r -^ J EtsyeMK^nASi) ' W 

As both sides of this equation are independent of the choice of the values of {o~i}, we 
may sum the right hand side over all possible values of {crj} and divide the result by 
2 N , where N is the total number of sites, to find 

1 



m 



2 N (2 cosh K„) N » 



J2J2 a 'i ex P( K P r ii a i°i) (Si)K- (29) 



PJ {ry}{<Ti} (ij) 



By inserting the identity 1 = Z(K P , {rij})/Z(K p , {%}) just in front of the sum over 
{(ij} in the summand, we obtain 

m = 2-(2coshA P )- E g g Wi) ^p^^) W* (30) 

= 2^2 cosh AF* ^ ^ exp(tf P 5>y*iO-i) (^(S^ (31) 
= (2 cosh i^* E «*p(*p 5>) fo)*^)*- (32) 



The last line fl32|) can be confirmed by applying the gauge transformation to ( |32| ) and 
using the same argument as above with gauge invariance of the product ((7j) K p (Si) k in 
mind to derive (0). Equation (^) is exactly the definition of the right hand side of 



fl), which completes the proof. A similar argument leads to (110) 
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